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MONTE  CARLO  CALCULATIONS  OF  LIGHT  SCATTERING  FROM  CLOUDS 


Gilbert  N.  Plass  and  George  W.  K&ttawar 
Abstract 

The  scattering  of  visible  light  by  clouds  is  calculated  from  an 
efficient  Monte  Carlo  code  which  follows  the  multiply  scattered  path 
of  the  photon.  The  single  scattering  phase  function  is  obtained  from  the 
Mie  theory  by  integration  over  a  particle  size  distribution.  The  photons 
are  followed  through  a  sufficient  number  of  collisions  and  reflections 
from  the  lower  surface  (which  may  have  any  desired  albedo)  until  they 
make  a  negligible  contribution  to  the  intensity.  Various  variance 
reduction  techniques  were  used  to  improve  the  statistics.  The  reflected 
and  transmitted  intensity  is  studied  as  a  function  of  solar  zenith  angle, 
optical  thickness,  and  surface  albedo.  The  downward  flux,  cloud  albedo, 
and  mean  optical  path  of  the  transmitted  and  reflected  photons  are  given 
as  a  function  of  these  same  parameters.  The  numerous  small  angle  scat¬ 
terings  of  the  photon  in  the  direction  of  the  incident  beam  are  followed 
accurately  and  produce  a  greater  penetration  into  the  cloud  than  is 
obtained  with  a  more  isotropic  and  less  realistic  phase  function. 
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INTRODUCTION 


The  multiple  scattering  of  a  photon  must  be  taken  into  account  in  any 
study  of  light  reflection  and  transmission  from  a  cloud.  Despite  the 
elegant  technique  developed  by  Chandrasekhar  [1960]  and  many  others,  it 
is  extremely  difficult  to  obtain  numerical  answers  when  the  single 
scattering  phase  function  is  strongly  anisotropic  and  when  the  properties 
of  the  atmosphere  change  appreciably  with  height.  Fritz  [1954,  1955] 
and  Tworoey  et  al  [1967]  have  obtained  interesting  results  pertaining 
to  clouds,  but  their  results  are  limited  to  larger  optical  depths  and 
directions  which  are  not  close  to  the  incident  direction  by  the  approxi¬ 
mations  inherent  in  their  methods. 

The  present  study  of  light  scattering  by  clouds  was  made  with  the 
Monte  Carlo  technique  which  has  now  become  practical  with  modern  high 
speed  digital  computers.  This  method  has  already  been  applied  to  many 
diverse  physical  and  mathematical  problems  [Hammersley  and  Handscomb.  1964], 
It  has  bten  used  by  Collins  and  Wells  [1965]  in  a  study  of  light  scattering 
by  atmospheric  aerosols. 

After  a  description  of  our  method  of  calculation,  the  reflected 
and  transmitted  light  from  a  cloud  is  discussed  as  a  function  of  optical 
thickness,  surface  albedo,  and  solar  incident  angle.  The  albedo  of  the 
cloud  is  presented  as  a  function  of  surface  albedo,  optical  thickness 
of  the  cloud,  and  incident  angle.  The  downward  flux  at  the  lower  surface 
of  the  cloud  and  the  mean  optical  path  of  both  the  reflected  and  trans¬ 
mitted  photons  is  discussed. 
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A  Lambert  reflecting  surface  and  a  single  distribution  of  sizes 
for  the  Mie  particles  is  assumed  in  all  these  calculations  in  order  to 
study  the  dependence  on  other  variables.  The  cloud  is  assumed  to  be 
horizontally  stratified  and  to  be  composed  of  water  droplets.  Effects 
due  to  incident  sky  radiation,  and  molecular  scattering  and  absorption 
art  not  considered  here.  The  num<  rous  small  angle  scatterings  of  the 
photon  in  the  direction  of  the  incident  beam  are  followed  accurately 
in  three  dimensions  in  the  present  calculation.  Among  the  interesting 
results  is  a  much  greater  penetration  of  the  beam  into  the  cloud  than 
is  obtained  with  a  more  isotropic  and  less  realistic  phase  function. 


MONTE  CARLO  METHOD 

The  Monte  Carlo  method  can  be  applied  to  any  problem  if  one  knows 
the  probability  for  each  step  in  a  sequence  of  events  and  desires  the 
probability  of  the  total  of  all  possible  events.  Thus  the  Monte  Carlo 
method  may  be  used  to  study  problems  in  radiative  transfer  [Collins  and 
Wells,  1965].  However,  some  of  the  methods  employed  in  this  work  diffei 
considerably  from  theirs. 

The  advantages  of  the  Monte  Carlo  method  for  the  study  of  radiative 
tranfer  in  planetary  atmospheres  are  as  follows: 

1)  The  calculations  may  be  performed  for  any  single  scattering 
phase  function  regardless  of  the  degree  of  anisotropy.  The  volume 
averaged  phase  function  may  he  extremely  anisotropic  for  aerosols 
and  water  droplets  in  the  atmosphere.  It  is  virtually  hopeless  to 
use  conventional  methods  relying  on  the  expansion  of  the  phase  function 
into  a  series  of  Legendre  polynomials  [Churchill  et  al,  1966;  Coulg  t 
etal,  1965], 


2)  Any  relevant  parameter  such  as  the  single  scattering  albedo, 
number  density  of  molecular  and  aerosol  scattering  centers,  and 
the  .ai xuus  cross-sections  may  vary  with  optical  depth  in  the 
atmosphere . 

3)  The  radiation  reflected  from  a  planetary  surface  may  follow 
any  desired  distribution,  instead  of  being  uniform  in  angle  as 
required  by  a  Lambert  surface.  The  light  reflected  from  many  real 
surfaces  varies  greatly  with  both  incident  and  reflected  angle 
[Coulson  et  al,  1965]. 

4)  A  single  computer  run  will  produce  fluxes  and  intensities  for 
both  a  large  number  of  surface  albedos  and  of  detectors  in  the 
atmosphere . 

5)  Any  reasonable  number  of  polar  and  azimuthal  angles  may  be 
selected. 

6)  The  average  mean  optical  path  of  both  the  reflected  and  transmitted 
photon  may  easily  be  calculated. 

The  disadvantage  of  the  Monte  Carlo  method  is  the  fact  that  the 
standard  deviations  of  the  results  are  roughly  inversely  proportional 
to  the  square  root  of  the  computing  times.  Therefore  the  method  may 
not  be  practical  when  high  accuracy  is  required. 

In  our  Monte  Carlo  code  the  three  dimensional  path  of  the  photon 
is  followed  through  a  cloud  and  the  atmosphere  as  it  is  multiply  scattered. 
The  single  scattering  phase  function  is  obtained  by  integration  of  the 
Mie  scattering  function  over  the  particle  size  distribution.  The  flux 
and  the  intensity  are  calculated  for  a  number  of  detectors  through  the 
atmosphere  for  a  number  of  surface  albedos.  In  the  present  calculations 
the  lower  surface  is  assumed  to  reflect  according  to  Lambert's  law. 
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Several  variance  reduction  techniques  minimized  the  computer  time. 
These  included  the  use  of  a  statistical  weight  for  the  photon  and  a 
requirement  for  forced  collisions  within  the  atmosphere.  A  new  method 
was  used  which  takes  account  of  all  orders  in  the  reflection  from  a 
Lambert's  surface  with  no  noticable  increase  in  the  computing  time. 

All  calculations  were  performed  on  an  IBM  360-50  computer.  The  random 
number  generator  was  the  multiplicative  congruential  generator  supplied 
with  the  scientific  subroutine  package.  When  using  this  generator  to 
choose  a  scattering  angle  from  the  cumulative  distribution  function  for 
Mie  scattering,  it  was  found  necessary  to  tabulate  the  function  at  very 
small  intervals  near  the  strong  forward  scattering  maximum.  A  linear 
interpolation  was  used  between  the  tabulated  values. 

Every  rejection  technique  used  in  the  program  was  tested  separately. 
However,  the  only  true  test  of  a  Monte  Carlo  program  is  in  the  results 
of  the  complete  program  itself.  The  reason  is  that  a  random  number 
generator  may  have  successive  random  numbers  that  have  the  right  properties, 
but  not  have  a  proper  distribution  of  n_  tuples.  Thus  our  results  were 
tested  in  a  number  of  cases  against  known  results  for  isotropic  and 
Rayleigh  scattering. 

The  Monte  Carlo  results  are  shown  in  Figure  1  together  with  the 
theoretical  results  of  Chandrasekhar  [1960]  for  isotropic  scattering 
from  a  semi-infinive  atmosphere  when  the  cosine  of  the  incident  angle 
uq  =  -0.2  and  the  single  scattering  albedo  u>o  =  0.2.  All  intensities  in 
this  paper  are  normalized  to  an  incident  intensity  of  unity  and  not  to  the 
value  it  as  is  sometimes  done.  All  fluxes  in  this  paper  are  normalized 


to  unit  incident  flux. 


The  Monte  Carlo  results  for  a  Rayleigh  phase  function  are  compared 
in  Figure  2  with  the  exact  solution  obtained  from  a  program  written  by 
Mr.  C.  N.  Adams  which  employs  invariant  imbedding  techniques.  In  this 
example  t  (the  optical  depth)  =  1  and  curves  are  shown  for  the  reflected 
and  transmitted  intensities  for  A  (surface  albedo)  =  0  and  A  =  0.8. 

The  agreement  in  all  cases  tested  appears  very  satisfactory. 


INTENSITY  REFLECTED  FROM  CLOUDS 

The  size  distribution  for  the  water  drops  in  a  cumulus  cloud  chosen 
for  this  calculation  is 

n(r )  =  2.373  r6  exp  (-1.5r)  ,  (1) 

where  the  radius  r  is  expressed  in  microns  and  the  concentration  n  is 
expressed  in  cm  u  .  This  is  the  distribution  function  used  by 
Deirmendjian  [1964],  We  assume  a  wavelength  of  0.7u  for  the  incident 
light  and  a  real  index  of  refraction  of  1.33  for  the  water  drops.  The 
results  are  also  applicable  to  other  visible  wavelengths  if  the  size 
of  the  particles  is  adjusted  in  proportion  to  the  wavelength. 

The  angular  intensity  function  for  single  scattering  is  calculated 
from  our  Mie  program  which  is  described  elsewhere  [Kattawar  and  Plass,  1967], 
The  results  are  averaged  over  the  size  distribution  of  the  drops  and  over 
the  two  directions  of  polarization.  The  resulting  intensity  function 
is  shown  in  Figure  3.  There  is  an  exceedingly  sharp  forward-scattering 
maximum  which  is  shown  in  more  detail  in  the  inset  in  the  upper  left  corner. 


From  this  function  a  cumulative  probability  distribution  was  calculated. 

It  was  necessary  to  specify  this  curve  at  over  240  points  with  variable 
spacing  in  order  to  be  certain  that  the  particles  are  scattered  according 
to  the  details  of  this  curve.  A  linear  interpolation  technique  was  used 
between  tabulated  points. 

The  reflected  and  transmitted  intensities  for  a  cloud  with  the  size 
distribution  given  by  Eq.  (1)  is  shown  in  Figures  4  and  5.  The  sunlight 
is  incident  vertically  so  that  uq  (cosine  of  incident  zenith  angle)  =  -1.0. 
Curves  are  given  for  clouds  with  an  optical  depth  t  from  0.01  to  30. 

The  intensity  is  shown  in  each  case  when  the  surface  albedo  A  =  0  and 
when  A  =  1.  The  transmitted  intensity  does  not  include  the  incident  photon 
until  it  has  been  scattered  for  the  first  time. 

When  a  cloud  is  sufficiently  thin,  the  scattered  intensity  approaches 
the  value  which  would  be  calculated  from  single  scattering  events  only. 

The  intensity  calculated  from  the  phase  function  shown  in  Figure  3  for 
single  scattering  only  is  indicated  in  Figures  4  and  5  by  small  squares 
when  t  =  0.01  and  0.1  and  when  A  =  0.  In  general  there  is  excellent 
agreement  between  these  results  and  those  of  the  Monte  Carlo  calculation. 

This  provides  another  check  on  the  Monte  Carlo  method.  Most  of  the  variation 
in  the  intensity  curve  as  a  function  of  u  is  due  to  the  nature  of  the 
phase  function  and  is  not  caused  by  fluctuations  In  the  result. 

The  single  scattering  va^ue  is  too  high  for  both  the  reflected  and 
transmitted  beam  for  the  range  0.0  <_  u  <_  0.1  and  t  =  0.1.  This  means 
that  the  cloud  is  already  sufficiently  thick  close  to  the  horizontal 
direction  when  t  =  0.1  so  that  multiple  scattering  must  be  taken  into  accouit. 
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For  the  range  0.9 <  u<_  1.0,  the  single  scattering  values  of  the  trans¬ 
mitted  intensity  are  too  low  for  both  x  =  0.1  and  0.01.  The  actual 
intensity  is  much  larger  because  of  the  many  small  angle  scattering 
events  which  occur  with  this  phase  function.  The  deviations  of  the 
multiply-scattered  intensities  from  the  single  scattering  values  for 
0.1  <_  u <_  0.9  provides  an  indication  of  the  statistical  fluctuation  in 
our  Monte  Carlo  results.  The  statistical  fluctuation  is  expected  to  be 
larger  when  the  intensities  are  small  (as  they  are  for  x  =  0.01)  rather 
than  when  they  are  closer  to  unity. 

When  the  surface  albedo  A  =  0,  the  reflected  intensity  as  shown 
by  Fig.  4  is  a  maximum  near  the  horizon  for  clouds  whose  optical  thickness 
x  is  small.  The  intensity  decreases  as  u""1  away  from  the  horizon  as  the 
number  of  single  scattering  centers  decreases  along  the  line  of  sight. 

The  maximum  in  the  curve  near  u  =  0.7  for  small  x  is  due  to  the  maximum 
in  the  angular  intensity  function  near  u  =  -0.75  as  shown  in  Figure  3. 

When  t  =  1,  the  reflected  intensity  is  more  nearly  uniform  at  all  angles 
because  of  multiple  scattering  effects.  For  large  x,  the  reflected 
intensity  has  a  minimum  at  the  horizon  and  a  maximum  at  the  zenith. 

When  A  =  1,  the  reflected  intensities  are  mucn  more  nearly  the 
same  at  all  angles.  The  intensity  is  essentially  uniform  over  angle 
for  small  x,  since  the  uniformly  reflected  radiation  from  the  Lambert 
surface  is  not  appreciably  modified  by  the  cloud.  For  clouds  with  large 
x,  the  reflected  intensity  is  a  minimum  near  the  horizon  since  the 
radiation  reflected  from  the  planetary  surface  can  not  penetrate  the 
cloud  at  these  large  zenith  angles. 


Tha  transmitted  intensity  in  Figure  5  has  a  large  max  mum  near  the 
zenith  for  thin  clouds  because  of  the  large  number  of  small  angle  scatterings 
around  the  original  direction  of  the  beam.  This  maximum  persists  even 
up  to  t  *  30  for  A  *  0,  although  it  is  very  much  less  pronounced.  There 
is  another  maximum  in  the  curves  for  thin  clouds  near  the  horizon  because 
of  the  cosine  effect.  However,  this  changes  to  a  minimum  for  thick  clouds. 

The  transmitted  intensity  as  shown  in  Figure  5  is  appreciably  larger 
near  the  horizon  at  all  t  for  A  =  1  than  for  A  =  0.  The  reason  is  that 
the  radiation  reflected  from  the  Lambert's  surface  near  the  horizon  has 
a  much  greater  chance  of  being  reflected  back  to  this  surface  if  it  makes 
one  or  more  small  angle  collisions  (see  Figure  3).  In  this  case  it  makes 
a  contribution  to  the  intensity  near  the  horizon.  In  order  to  contribute 
to  the  downward  intensity  at  an  angle  close  to  the  zenith,  the  photon 
must  undergo  large  angle  scattering  which  is  much  less  probable. 

When  =  1.0,  the  intensities  are  azimuth  independent.  Thi ;  is 
not  the  case  for  other  angles  of  incioence.  We  can  use  as  many  intervals 
as  desired  for  both  the  zenith  and  azimuthal  angles  in  our  Monte  Carlo 
program  within  certain  restrictions  on  memory  size.  It  was  decided  to 
use  only  two  intervals  for  the  azimuthal  angle  for  other  angles  of 
incidence  for  ease  of  presentation  of  the  results.  The  results  for 
uq  *  -0.5  are  shown  in  Figs.  6-9.  Take  $  =  0°  at  the  intersection  of 
the  horizontal  plane  and  the  plane  containing  the  incident  beam  and 
in  the  direction  of  the  incident  beam.  Then  the  #  intervals  used  in  this 
calculation  are:  (1)  -:c°  <_  *  <_  90° ;  (2)  90°  <_  <*>  <_  270°.  The  results 
for  interval  (1)  are  shown  on  the  left  hand  side  of  Figures  6-°  and 
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those  for  interval  (2)  on  the  right  hand  side,  i.e.  one  intensity  curve 
shows  the  variation  from  one  horizon  to  the  zenith  back  to  the  other 
horizon  averaged  in  each  case  over  the  4  angles  indicated  above. 

The  reflected  intensity  for  thin  clouds  shown  in  Figure  6  for 
A  =  0  is  largest  near  the  horizon  nearest  the  direction  of  the  incident 
beam.  It  falls  off  very  regularly  to  a  minimum  near  the  zenith  and  then 
rises  less  rapidly  to  a  smaller  value  on  the  opposite  horizon.  The  value 
at  the  minimum  is  500  times  smaller  than  at  the  maximum.  This  curve 
results  both  from  the  cosine  term  near  the  horizon  and  the  greatly  reduced 
probability  for  backscattering  compared  to  forward  scattering.  As  the 
optical  depth  of  the  cloud  increases,  the  intensity  curves  follow  qualita¬ 
tively  the  same  variation,  but  the  magnitude  of  the  variations  is  greatly 
reduced  by  multiple  scattering. 

The  intensity  when  A  =  1  again  has  relatively  little  dependence 
on  angle  as  Is  shown  in  Figure  7.  The  intensity  is  nearly  constant  for 
thin  clouds.  For  thick  cxouds  it  has  a  small  maximum  near  the  point 
on  the  horizon  closest  to  the  direction  of  the  original  beam  and  a  minimum 
at  the  opposite  horizon. 

When  A  =  0  the  transmitted  intensity  as  shown  in  Figure  8  has  a 
maximum  in  the  direction  of  the  original  beam  for  t  <_  3.  This  is  expected 
from  the  numerous  multiple  small  angle  scattering  events  around  the 
original  direction.  The  minimum  transmitted  intensity  is  near  0  -  0.8 
on  the  opposite  side  from  the  direction  of  the  original  beam.  This  is 
caused  by  the  single  scattering  probability  decreasing  rapidly  as  the  angle 
increases  in  this  range.  For  thick  clouds  the  intensity  is  relatively  constant, 
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but  shows  a  tendency  to  be  two  to  three  times  smaller  near  the  horizon 
than  at  the  zenith. 

When  A  =  1,  the  qualitative  shape  of  the  curves  shown  in  Figure  9 
is  the  same  as  in  Figure  8.  The  values  are  larger  near  the  horizons 
as  the  albedo  increases  due  to  the  photons  reflected  from  the  Lambert 
surface  into  an  angle  near  the  horizon  and  then  back-scattered  into  a 
downward  direction  near  the  horizon  by  the  many  probable  small  angle 
scatterings.  The  curves  for  large  x  are  very  uniform. 

Similar  curves  for  sunlight  incident  near  the  horizontal  at  uq  =  -0.1 

are  given  in  Figures  10-13.  The  reflected  intensity  shown  in  Figure  10 

is  similar  to  the  curves  for  u  =  -0.5  except  that  the  variations  in 

o 

the  curves  are  greater  and  all  curves  show  a  pronounced  minimum  near 
the  zenith.  Even  a  thick  cloud  reflects  ten  times  as  much  light  toward 
the  horizon  in  the  direction  of  the  incident  beam  as  it  does  to  the 
zenith.  The  photons  can  escape  from  the  cloud  easily  in  this  direction 
by  many  probable  small  angle  scatterings. 

The  reflected  intensity  when  A  =  1  also  shows  more  variation  than 
before  as  shown  in  Figure  11.  All  curves  have  a  maximum  at  the  near 
horizon  and  a  minimum  near  the  zenith. 

The  transmitted  intensity  as  shown  in  Figure  12  for  the  case  A  -  0 
shows  the  expected  maximum  for  thin  cJouds  around  the  direction  of  the 
original  beam.  For  thin  clouds  the  maximun  transmitted  intensity  at  the 
horizon  is  6000  times  larger  than  the  minimum  value  near  the  zenitn. 

When  A  =  1,  Figure  13  shows  again  that  the  effect  of  reflected  photons 
from  the  planetary  surface  is  predominately  to  increase  the  transmitted 
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intensity  near  the  horizon.  It  is  interesting  to  note  in  most  cases 
when  t  <_  1  that  there  is  a  considerable  region  in  which  the  transmitted 
intensity  is  nearly  proportional  to  exp  [-(const )y].  When  t  is  large 
the  intensity  is  several  times  larger  at  the  zenith  than  at  the  horizon. 

The  variation  of  the  reflected  and  transmitted  intensities  with 
solar  angle  is  further  illustrated  by  figures  14  and  15  for  =  -0.02. 

When  the  incident  beam  is  so  close  to  the  horizon,  the  reflected  intensity 
has  a  very  strong  maximum  at  the  forward  horizon  even  for  thick  clouds. 

When  t  =  10,  the  intensity  near  the  horizon  is  more  than  100  times  greater 
than  that  at  the  zenith  and  20  times  greater  than  the  intensity  at  the 
far  horizon.  The  shapes  of  the  transmitted  intensity  curves  are  rather 
similar  for  uq  =  -0.02  and  =  -0.1.  However,  the  transmitted  intensity 
is  smaller  in  the  former  case  because  of  the  larger  cloud  albedo. 

Some  of  the  results  of  Twomey  et  al  [1067]  may  be  compared  with  these 
curves.  This  method  is  only  valid  for  i  >  1  and  does  not  take  account  '"f 
the  numerous  small  angle  scatterings  accurately  as  is  done  in  our  calculation. 
Our  curves  show  in  general  stronger  maxima  around  the  direction  of  the 
original  beam  and  differ  numerically  from  some  of  his  results. 


rLux 

The  downward  diffuse  flux  at  the  lower  boundary  of  the  cloud  as  a 
function  of  optical  thickness  is  shown  in  Figure  16  for  A  =  0.  A  photon 
is  not  counted  in  the  difruse  flux  until  it  lias  undergone  one  scattering 
event.  All  fluxes  are  normalized  to  unit  incident  flux.  The  flux  is 
small  for  small  t,  since  there  are  not  enough  water  droplets  to  scatter 


an  appreciable  number  of  photons.  The  flux  increases  with  t  until  the 
curves  reach  a  maximum.  The  intensity  decreases  rapidly  with  further 
increases  in  t  ,  since  the  cloud  is  so  thick  that  the  photons  have  a 
greater  probability  of  escaping  from  the  upper  surface  tf  an  in  finding 
their  way  to  the  bottom.  The  maximum  flux  occurs  at  successively  smaller 
values  as  the  incident  beam  approachs  the  horizon. 

When  the  surface  albedo  is  high,  the  downward  flux  at  the  lower 
surface  of  the  cloud  under  certain  conditions  may  be  greater  than  the 
incident  flux,  as  was  first  pointed  out  by  fritz  [ 19 5 b J .  The  downward 
flux  (including  the  contribution  from  the  unscattered  incident  flux 
normalized  to  unity  at  the  cloud  top)  is  shown  in  Figure  17  as  a  function 
of  optical  thickness  tor  four  different  values  of  when  A  =  1.  When 

U  =  -1.0  and  c  =  6,  the  downward  flux  is  22'L  greater  than  the  incident 

o 

flux.  When  u  =  -0.5,  the  downward  flux  is  0.13%  greater  than  the  incident 
o 

flux  when  t  =  0.1.  The  downward  flux  is  always  less  than  unity  for 
near  the  horizon.  The  large  effect  wh*“:i  the  incident  beam  is  vertical 
is  connected  with  the  strong  forward  scattering  of  the  particles.  A 
relatively  large  fraction  of  the  radiation  is  transmitted  through  the 
cloud  near  the  vertical,  but  finds  it  difficult  to  escape  back  out  because 
the  photons  are  now  traveling  in  almost  random  directions.  When  UQ  -  -1.0 
and  T  =  0,  the  downward  fJux  is  greater'  than  unity  for  all  A  >  0.4.  Tor 
example,  the  downward  flux  is  1.0994  when  A  =  1;  1.0bb2  when  A  =  0.8; 

1.0350  when  A  =  O.b;  0.99^7  when  A  -  0.4;  0.9514  when  A  =  0. 


The  strong  penetration  of  radiation  into  a  cloud  composed  of  water 
droplets  larger  than  the  wavelength  is  shown  by  a  comparison  of  the  flux 
calculated  here  with  that  for  Rayleigh  scattering.  For  example,  when 
=  -1.0  and  A  =  1  the  downward  flux  is  0.6573  for  Rayleigh  scattering, 
but  1.0994  for  Mie  scattering.  The  difference  is  even  larger  between 
the  diffuse  intensity  in  the  zenith  direction  which  is  0.0977  for  Rayleigh 
scattering  and  0.648  for  Mie  scattering. 


MEAN  0PTTCAL  PATH 

The  mean  optical  path  of  the  photon  is  a  quantity  of  interest  in 
several  types  of  problems.  For  example,  when  the  reflected  light  from 
clouds  is  measured  from  a  satellite,  a  knowledge  of  the  mean  optical  path 
is  necessary  to  allow  for  molecular  absorption  at  certain  wavelengths. 

The  mean  optical  path  of  the  reflected  photon  is  shown  in  Figure  18 

for  various  angles  of  incidence.  For  thin  clouds  of  a  given  optical 

thickness,  the  mean  optical  path  is  greater  the  closer  the  incident  beam 

is  to  the  horizon.  The  pi  oton  in  this  case  travels  a  greater  distance 

along  the  beau  direction  before  being  scattered  back  out. 

For  thick  clouds  of  a  given  optical  thickness,  the  mean  optical  path 

is  greatest  for  an  incident  beam  at  the  zenith.  In  this  case  the  photon 

on  the  average  penetrates  deeply  into  the  cloud.  Because  of  the  single 

scattering  phase  function,  the  photon  is  more  likely  to  be  scattered 

in  an  upward  direction  by  a  single  collision  when  the  incident  beam  is 

nearly  horizontal.  For  example,  when  the  cloud's  optical  thickness 

t  =  10,  the  mean  optical  path  for  the  reflected  photon  is  24.2  for 

u  =  -1.0;  20.8  for  u  =  -0.5;  12.9  for  u  =  -0.1;  9.47  for  u  =  -0.02. 
o  o  o  o 
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The  curves  cross  over  one  another  in  the  region  from  t  =  1  to  10. 

The  mean  optical  path  of  the  transmitted  photon  is  shown  in 
Figure  19.  This  quantity  for  a  given  cloud  thickness  always  increases 
as  the  incident  beam  moves  from  the  zenith  to  the  horizon.  The  differences 
for  various  incident  angles  are  large  for  thin  clouds  and  very  small 
for  thick  clouds.  When  uq  =  -1.0,  the  photon  efficiently  penetrates 
the  cloud  as  it  undergoes  many  small  angle  forward  scatterings.  When 
Uq  a  -0.02,  the  photon  on  the  average  travels  a  distance  of  the  order 
t  *  1  before  it  undergoes  a  large  angle  scattering  collision.  After  one 
large  angle  collision  the  photon  essentially  loses  its  memory  of  its 
initial  direction.  Thus  all  photons  travel  approximately  the  same  distance 
from  this  point  until  they  leave  the  lower  surface  of  the  cloud.  This 
is  the  reason  all  the  mean  optical  paths  are  nearly  the  same  for  large 
t.  For  example,  when  t  =  10,  the  mean  optical  thickness  of  the  trans¬ 
mitted  photon  is  19.2  for  u  -  -1.0;  24.9  for  u  =  -0.5;  28.0  for  u  =  -0.1, 

o  o  o 

28.2  for  u  =  -0.02. 
o 

Mean  optical  paths  for  t  >  2  and  uq  =  -0.55  have  been  given  by 
Twomey  et  al  [1967],  The  curves  are  qualitatively  similar  to  those 
shown  here.  Our  reflected  mean  rptical  paths  are  larger  than  theirs. 

This  is  presumably  because  our  calculations  accurately  allows  for  the 
numerous  small  angle  scatterings  which  allow  the  photon  to  penetrate 
deeper  into  the  cloud  and  thus  to  have  a  greater  optical  path.  It 
may  also  be  due  to  the  fact  that  their  single  scattering  phase  function 


is  somewhat  different  from  ours. 
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CLOUD  ALBEDO 


The  cloud  albedo  is  defined  as  the  fraction  of  the  incident  light 
reflected  from  the  cloud.  It  is  shown  in  Tigures  20-24  for  various  angles 

of  incidence  u  and  surface  albedos  A. 

o 

The  cloud  albedo  for  =  -1.0  is  shown  in  Figure  20.  When  the 
planetary  surface  is  completely  absorbing  (A  =  0),  the  cloud  does  not 
reflect  appreciable  radiation  until  i  is  of  the  order  of  unity.  The 
cloud  albedo  then  rises  rapidly  until  it  is  0.88  when  t  =  100.  When 
A  =  1,  the  cloud  albedo  has  the  value  unity  for  all  t,  since  no  absorption 
was  assumed  by  the  water  droplets  in  this  calculation. 

For  A  =  0.8,  the  cloud  albedo  is  0.8  when  the  cloud  is  sufficiently 
thin.  As  t  increases,  the  albedo  first  decreases  since  the  increased 
reflection  by  the  cloud  is  more  than  offset  by  the  increased  downward 
scattering  of  the  radiation  reflected  form  the  planetary  surface.  The 
minimum  cloud  albedo  in  ti.is  case  is  0.782  near  t  '•  3.  A  similar  minimum 
also  occurs  for  A  =  0.1*  and  0.6.  The  cloud  albedo  rises  as  t  inci  eases 
beyond  this  point. 

The  cloud  albedos  for  u  =  -O.b,  -0.1,  and  -0.02  are  shown  in 

o 

Figures  21-23.  These  curves  have  the  same  shape  as  that  for  uq  =  -1.0, 

except  that  there  is  no  longer  a  small  minimum  in  t lie  curves  for  A  =  0.4 

to  0.8.  Because  of  the  strong  forward  scattering,  the  cloud  albeuo  always 

increases  with  t  when  I  u  I  <  0.5. 

'o'  — 

The  cloud  albedo  also  always  increases  as  |n^|  decreases  for  a 

fixed  t  and  A.  For  example,  when  t  =  0.1  and  A  =  0,  the  cloud  albedo 

is  0.00468  when  u  =-1.0;  0.0181  when  u  =  -O.b,  0.180  when  u  =  -0.1; 

o  o  o 

|j 

)' 


0.492  when  pQ  s  -0.02.  When  t  *  10  end  A  =  0,  the  cloud  albedo  increases 

from  0.466  when  u  =  -1.0  to  0.852  when  p  =  -0.02. 

o  o 

The  cloud  albedos  shown  here  are  appreciably  smaller  than  those 
given  by  Fritz  [1954]  in  his  Figure  4.  This  is  probably  explained  by 
the  penetration  to  a  greater  depth  by  the  actual  photons  through  multiple 
small  angle  scattering  as  described  in  our  calculation.  Our  values  are 
also  usually  less  than  those  of  Korb  and  Holler  [1962]  presumably  for  the 
same  reason. 

The  curves  of  cloud  albedo  for  various  values  of  A  approach  each 

other  more  closely  as  t  increases.  This  merely  states  that  a  relatively 

smaller  percentage  of  the  photons  is  able  to  reach  the  planetary  surface 

and  return  to  the  upper  surface  as  t  becomes  large.  An  important  point 

is  that  there  can  be  significant  variations  in  the  cloud  albedo  as  the 

surface  albedo  changes  for  values  of  x  between  10  and  100,  as  shown  in 

Figures  20-23.  The  optical  thickness  of  Venus  is  usually  placed  in  this 

range.  Table  I  shows  the  cloud  albedo  for  t  =  10  and  30,  for  various 

values  of  u  and  surface  albedo  A.  It  is  immediitely  obvious  that  the 
o 

reflected  flux  still  depends  appreciably  on  the  surface  albedo  even  for 
large  t.  Thus  there  is  the  possibility  of  obtaining  information  about 
conditions  beneath  a  thick  cloud  layer  from  a  study  of  the  reflected 
light  because  of  the  strong  penetration  in  the  direction  of  the  original 
beam  allowed  by  the  many  small  angle  scatterings  undergone  by  the  photon. 


TABLE  I.  CLOUD  ALBEDO 


T 

Uo 

A  =  0 

A  =  0 

10 

-1.0 

0.466 

0.519 

10 

-0.5 

.63d 

.674 

10 

-0.1 

.799 

.819 

10 

-0.02 

.852 

.867 

30 

-1.0 

.773 

.786 

30 

-O.b 

.818 

.827 

30 

-0.1 

.900 

.905 

A  =  0.4 

A  =  0.6 

o 

II 

< 

0.588 

0.679 

0.807 

.721 

.782 

.869 

.845 

.880 

.928 

.886 

.912 

.947 

.803 

.831 

.881 

.840 

.861 

.901 

.912 

.924 

.946 
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LEGENDS  FOR  FIGURES 

Reflected  intensity  from  a  semi-infinite  atmosphere  with  isotropic 
scattering,  mq  (cosine  of  incident  angle)  =  -0.2,  and 
(single  scattering  albedo)  -  0.2  as  a  function  of 
u  (cosine  of  zenith  angle).  The  continuous  curve  is  calculated 
from  Chandrasekhar  [I960].  The  Monte  Carlo  results  are  for 
20,000  histories.  The  incident  intensity  is  normalised  to  unity. 
Reflected  and  transmitted  intensity  for  a  Rayleigh  phase 
function  for  an  atmosphere  with  t  (optical  depth)  =  1.  The 
upper  curves  are  for  A  (surface  albedo)  =  0.8  and  the  lower 
curves  are  for  A  =  0.  The  continuous  curve  was  calculated 
by  a  program  written  by  Mr.  C.  N.  Adams.  The  Monte  Carlo 
results  are  for  30,000  histories. 

Angular  intensity  function  for  Mie  scattering  as  a  function  of 
the  cosine  u  of  scattering  angle  averaged  over  the  size  distri¬ 
bution  given  by  Eq.  (1)  and  over  the  two  directions  of  polarization. 
The  inset  in  upper  left  shows  the  curve  near  u  -  1.  It  is  assumed 
that  the  wavelength  of  the  incident  light  is  0.7y  and  that  the 
index  of  refraction  of  the  water  drops  is  1.33. 

The  reflected  intensity  is  shown  as  a  function  of  u,  the  cosine 
of  the  zenith  angle.  The  curves  on  the  left  and  right  portion 
of  the  figure  are  for  A  (surface  albedo)  =  0  and  1  respectively. 

The  optical  depth  of  the  cloud  is  t.  The  sunlight  is  incident 
vertically,  uq  (cosine  of  incident  zenith  angle)  =  -1.0.  The 
incident  intensity  is  normalized  to  unity. 
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Fig.  5.  Transmitted  intensity  for  =  -1.0  and  A  =  0  and  1.  The 
transmitted  intensity  does  not  include  the  incident  photon 
until  it  has  been  scattered  the  first  time.  S'.e  caption 
for  Fig.  4. 

Fig.  6.  Reflected  intensity  for  Uq  =  -0.5  and  A  =  0  as  a  function 

of  y,  the  cosine  of  the  zenith  angle.  The  left  hand  portion 
of  the  graph  refers  to  values  averaged  over  the  azimuthal  angle 
for  90°  on  both  sides  of  the  direction  of  the  original  team. 

The  values  on  the  right  portion  of  the  graph  are  for  values 
averaged  ov»r  the  remaining  azimuthal  angles.  Thus  one 
intensity  curve  from  left  to  right  she  s  the  variation  from 
one  horizon  to  th«  zenith  and  back  to  the  other  horizon  averaged 
over  the  indicated  azimuthal  angles. 

Fig.  7.  Reflected  intensity  for  Uq  =  -0.5  and  A  =  1.  See  caption  for  Fig.  6 

Fig.  8.  Transmitted  intensity  for  =  -0.5  and  A  *  0.  See  caption 

for  Fig.  6 . 

Fig.  9.  Transmitted  intensity  for  =  -0.5  and  A  =  1.  See  caption 

for  Fig.  6. 

Fig-  10.  Reflected  intensity  for  u0  =  -0.1  and  A  =  0.  See  caption  for  Fig.  b 

Fig.  11.  Reflected  intensity  for  =  -0.1  and  A  =  1.  See  caption  tor  Fig.  6 

Fig.  12.  Transmitted  intensity  for  =  -0.1  and  A  =  0.  See  caption 

for  Fig.  6. 

Fig.  13.  Transmitted  intensity  for  =  -0.1  and  A  =  1.  See  caption 

tut  Fig.  6. 

Fig.  14.  Reflected  intensity  for  =  -0.02  and  A  =  0.  See  caption 


for  Fig.  6. 
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Fig.  15.  Transmitted  intensity  for  uq  =  -0.02  and  A  =  0.  See  caption 
for  Fig.  6. 

I  i 

Fig.  16.  Downward  diffuse  flux  at  lower  boundary  of  cloud  as  a  function 

of  optical  thickness  for  uq  =  -1.0,  -0.5,  -0.1  and  for  A  =  0. 

The  incident  flux  is  normalized  to  unity  in  each  case. 

Fig.  17.  Downward  diffuse  flux  plus  remaining  flux  from  incident  beam 

at  lower  boundary  of  cloud  as  a  function  of  optical  thickness 

of  cloud  for  u  =  -1.0,  -0.5,  -0.1,  -0.02  and  for  A  =  1.  The 
o 

incident  flux  is  normalized  to  unity  in  each  case. 

Fig.  18.  Mean  optical  path  of  reflected  photon  as  a  function  of  optical 

thickness  of  cloud  for  =  -0.02,  -0.1,  -0.5,  -1.0. 

Fig.  19.  Mean  optical  path  of  transmitted  photon  as  a  function  of  optical 

thickness  of  cloud  for  u  =  -0.02,  -0.1,  -0.5,  -1.0. 

o 

Fig.  20.  Cloud  albedo  as  a  function  of  optical  thickness  for  uq  =  -1.0 
and  A  =  0,  0.2,  0.4,  O.b,  0.8. 

Fig.  21.  Cloud  albedo  as  a  function  of  optical  thickness  for  mq  =  -0.5 
and  A  =  0,  0.2,  0.4,  0.6,  0.8. 

Fig-  22.  Cloud  albedo  as  a  function  of  optical  thickness  for  Uq  =  -0.1 
and  A  =  0,  0.2,  0.4,  0.6,  0.8. 

Fig.  23.  Cloud  albedo  as  a  function  of  optical  thickness  for  =  -0.02 
and  A  =  0,  0.2,  0.4,  0.6,  0.8. 
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The  scattering  of  visible  light  by  clouds  is  calculated  from  an  efficient 
Monte  Carlo  code  which  fellows  the  multiply  scattered  path  of  the  photon.  The 
single  scattering  phase  function  is  obtained  from  the  Mie  theory  by  integration 
over  a  particle  size  distribution.  The  photons  are  followed  through  a  sufficient 
number  of  collisions  and  reflections  from  the  lower  surface  (which  may  have  any 
desired  albedo)  until  they  make  a  negligiblt  ventribution  to  the  intensity. 

Various  variance  reduction  techniques  were  used  to  improve  the  statistics.  The 
reflected  and  transmitted  intensity  is  studied  as  a  function  of  solar  zenith  angle, 
optical  thickness,  and  .surface  albedo.  The  downward  flux,  cloud  albedo,  and  mean 
optical  path  of  the  transmitted  and  reflected  phe tons  are  given  as  a  function 
of  these  same  parameters.  The  numerous  small  angle  scatterings  of  the  photon  in 
the  direction  of  the  incident  beam  are  followed  accurately  and  produce  a  greater 
penetration  in.o  the  cloud  than  is  obtained  with  a  more  isotropic  and  less 
realistic  phase  function. 
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